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From counts of 60000 stars in the neighborhood of the Fornax dwarf galaxy, the distribution of the 
projected stellar density in the galaxy has been derived. It is found to be smooth, with its ellipticity aver- 
aging about 0.35, and with some asymmetry in the central region. The density profile falls off more rapidly 
with distance than Hubble’s interpolation formula for giant ellipticals. This is interpreted as the result 
of the gravitational influence of the Galaxy on the outer stars of the Fornax*dwarf. From estimates of the 
mass of Fornax it is found that the predicted tidal limiting radius is in reasonable agreement with that 
observed. There does not appear to be any general obscuration within the galaxy; distant galaxies are 


distributed uniformly behind it. 


I, INTRODUCTION 


HE distribution of luminosity with radius for 
elliptical galaxies has been the subject of numer- 
jus studies (Hubble 1930; Redman and Shirley 1938; 
e Vaucouleurs 1948; foam 1956; Liller 1960; sa 
any others). In general, these investigators have 
ound that the interpolation formula employed by 
ubble fits the intensity distribution of giant elliptical 
alaxies well, although there seems to be a tendency for 
he luminosity to fall off more rapidly than the law at 
‘reat distances. Almost without exception the elliptical 
ralaxies so far examined photometrically, of which there 
re some 50, have been giant objects of high luminosity 
ith bright nuclei and steep luminosity gradients. Two 
xceptions are NGC 205 and M32, which have been 
measured by Redman and Shirley (1938), Sytinskaya 
1943), and van Houten, Oort, and Hiltner (1954). 
3oth of these galaxies have moderate absolute mag- 
itudes, bright nuclei, and fairly steep luminosity 
rradients. There are no published measures, however, 
of the kind of elliptical that appears to be the most 
bundant type, that is, the Sculptor-type dwarf 
sllipticals, which have faint absolute magnitudes, no 
uclei, and very low luminosity gradients. 

This paper reports the examination of the structure 
of one such galaxy, the Fornax dwarf [a=2'38™, 
j= —34°7 (1950)], discovered by Shapley (1938) 
Fig. 1). It is apparently the largest and most populous 
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of the Sculptor-type galaxies in the local group, and 
the only one to contain globular clusters (Baade and 
Hubble 1939; Shapley 1939; Hodge 1961). Its distance 
is approximately 188 kpc, as determined by Baade and 
Hubble (1939), who measured the brightest stars in 
two of the globular clusters by photographic transfers. 
(There are three sources of error in this determination: 
the errors in the primary sequence used, the errors due 
to the photographic transfer from a sequence to the very 
low altitude reached by Fornax at Mount Wilson, and 
the uncertainty in the identification of the magnitude of 
the brightest stars in Fornax’s globular clusters with 
any particular absolute magnitudes. These total at 
least 0.5 magnitude and probably more.) 

The Fornax dwarf galaxy has such an extremely low 
surface brightness that conventional methods of deter- 
mining the distribution of intensity are relatively ineffec- 
tive. Photoelectric scans across the face of the Fornax 
dwarf have not yet been successful. Small-camera photo- 
graphs can be obtained which show the galaxy, but only 
the central portion. Large reflectors can detect the indi- 
vidual brightest stars, but the small field limits them 
to only a fraction of the galaxy on any one plate. 
Therefore, the most effective method of determining 
the structure of such an object with present equipment 
is by means of star counts on Schmidt plates, which 
have large enough fields to show the whole galaxy and 
a large area around it, and which are capable of very 
faint limiting magnitudes. 
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Fic. 1. The Fornax dwarf galaxy and its five globular clusters (see paper I). Photographed in the blue by the ADH Baker- Schma 


The major axis runs nearly VE by SW, and extends almost to the corners of the photograph. 


Il. OBSERVATIONAL DATA 


Two sets of Schmidt plates have been taken of the 
Fornax dwarf for purposes of making star counts: one 
set with the ADH Schmidt at the Boyden Observatory, 
South Africa, and one with the 48-inch Schmidt of the 
Palomar Observatory. The low altitude reached by 
Fornax at Palomar was found to offset the advantage 
of the superior aperture of the 48-inch telescope, so 
that the ADH plates on the best nights had the smaller 
images and the fainter limiting magnitudes. This report 


TaBLE I. Observational material. 


Total 
Plate Exposure stars Background 
number Aperture (minutes) counted (stars per []’) 
ADH 71 33-inch 75 26 600 0.70 
ADH 5101 24-inch 95 12 700 0.56 
_ ADH 5094 24-inch 90 11 100 0.52 
ADH 5089 24-inch 60 9650 0.50 


is based on counts made on four ADH plates take 
under superior seeing conditions (Table I). Three o| 
the plates were exposed with the 24-inch diaphragn 
in place, so that the vignetting-free field for these plate; 
is 6 inches in diameter. The remaining plate, ADH 71 
was exposed without the diaphragm (aperture: 3< 
inches) in order to reach a fainter limiting magnitude 
For this plate the vignetting-free field is only 2 inche: 
in diameter. Tests made with the 24-inch aperture in 
place have shown that there is no systematic distanc¢ 
correction within the 6-inch diameter vignetting-fre¢ 
field (Hodge 1960). For plates at full aperture, the 
field corrections have been determined by Lindsay 
(1953). They are large enough that plate ADH 71 car 
be used only as an approximate indicator for the outel 
star distribution. 

For the star counts a 2.5-mm grid was superimposed 
over the plate (scale: 68’ per mm), oriented VS and 
EW. The number of stars was tabulated for each square 
over an area extending 80 mm in each direction, with 
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Fic. 2. Equidensity contours for plate ADH 71. Crossed circles 
are globular clusters. Density is in unit of stars per counting 
square (8 square minutes of arc). The outer contour line represents 
a value of 10 stars per square, and the density increases inwards 
in steps of 5 stars per square. 


several auxiliary areas extending 40 mm further. This 
provided a complete count for an area of some 12 000 
square millimeters, or 15000 square minutes of arc 
(4.2 square degrees). Plate ADH 71 was counted three 
times, the remainder each counted once. Counts were 
made to the limiting magnitude, and each complete 
count was made at one long sitting. The latter was 
found necessary in order to avoid systematic effects. 
The three counts of ADH 71 agree well, and a 5% 
mean error represents the reliability of the count of any 
one square. The counts can be considered internally 
quite consistent. 

Only the long-exposure plate, ADH 71, had images 
crowded together enough to make counting difficult. 
For the few squares at the center of the galaxy, the 
images on this plate were quite crowded, and the final 
counts uncertain. The plate ADH 5101 provided similar 
total counts (some 50 stars per square) at the center, 
but the images were smaller than for ADH 71, so that 
no difficulty was experienced, and the central counts 
for this plate can be considered reliable. 


Fic. 3. Contours for plate ADH 5101. Outer line has 
same density as for ADH 71. 
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Fic. 4. Contours for plate ADH 5094. Outer line has 
same density as for ADH 71. 


III, DISTRIBUTION OF STARS 
Lines of Equidensity 


In order to represent the counts graphically, equi- 
density contours have been constructed from the data 
for each plate. This was done by plotting positions of 
equal counts and drawing a smooth curve through all 
the points for (v+2)+2, with n=10, 15, 20---, where x 
is the number of stars counted per square. For example, 
the equidensity line for 7=15 is as plotted actually the 
line for the average position of all squares with »=15 
through n=19. The equidensity contour diagrams 
derived in this way from the four plates are shown in 
Figs. 2 through 5. These differ from each other some- 
what in detail, but show the same general structure for 
the galaxy. The differences are due to the different 
limiting magnitudes, the certain amount of arbitrariness 
involved in the construction of the outermost and the 
innermost contours and any small irregularities in 
emulsion sensitivity. (The very low density gradient 
of the galaxy and the small steps between equidensity 
contours tend to exaggerate statistical fluctuations and 


Fic. 5. Contours for plate ADH 5089. Outer line has 
same density as for ADH 71. 


DSW PUGH IL. 
0.4 
€ cons 
One 
Co) 10 20 30 40 
r 


Fic. 6. Variation of the ellipticity of the Fornax galaxy with 
distance from the center (7 is in units of minutes of arc). 


irregularities far more than for isophototal maps of 
normal elliptical glaxies.) The counts and contours show 
that there are no real irregularities in the galaxy; the 
distribution of stars is smooth. 


Orientation 


The equidensity contour diagrams provide the best 
means of determining the orientation of the galaxy. 
The value obtained by averaging the orientations for 
all contours in Figs. 2 through 5 is 54° for the major 
axis, measured eastward from the north. There is no 
systematic change of orientation with radius, though 
some of the drawn contours might indicate such a 
situation. 

Ellipticity 


As in giant elliptical galaxies (Liller 1960), the ellip- 
ticity of the Fornax galaxy contours varies with the 
distance, increasing outwards. The ellipticity is defined 
here by the standard relation (Liller 1960) 


e=1—(0/a), (1) 


where a is the major and 6 the minor axis. Figure 6 
shows how the ellipticity varies with distance along the 
major axis. Because of the asymmetry of the galaxy, 
however, the ellipticities derived are not entirely 


TaBLeE II. Density profiles—major and minor axes.* 
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realistic, and must be treated as approximate indicator 
of the true situation, which cannot accurately be repre 
sented by elliptical parameters. 


Density Profiles 


Stellar densities along the major and minor axes ar} 
given in Table II. The values are averages for fou) 
counting squares, reduced to the number of stars pe 
square minute. The area represented by each entry 1 
Table II is 32 sq minutes, so that the statistical uncer 
tainty of each value of the density of stars per squar 
minute, is ¥ p/32. Major-axis density profiles are show} 
in Figs. 7 and 8. Here it is evident that on the long 
exposure plate, ADH 71, saturation of images wa 
reached in the center.of the Fornax galaxy. 

For purposes of determining the distribution of star} 
in the outermost portions of Fornax, and in thi 
foreground, east-west profiles to greater distances fron 
the center of the galaxy were obtained. The counts ar 
shown in Fig. 9, along with mean curves fitted by eye 
Adopted values for the density of foreground stars oi 
each plate are given in Table I. Apparently becaus| 
of the high galactic latitude of the Fornax dwar) 
(b= —64°), there was no detectable variation of th} 
foreground count across that portion of the plate wher} 
counts were made. 


Asymmetry 


Figures 7 and 8 show that there is a striking asym} 
metry in the major-axis profiles, evident from all fou} 
plates counted. This asymmetry has the characteristic} 
of a displacement of the central, highest density poin} 
by about 6’ from the geometrical center of the galaxy } 
determined by the outer regions. 

It is difficult to explain this asymmetry in a straight] 
forward way. If the age of the galaxy is less than, or oj 
the order of, the relaxation time, then the presen} 
asymmetry might. reflect the situation in the early) 
stages of its development. The age of the galaxy 1 


ADH 71 ADH 5101 ADH 5094 ADH 5089 iM 

r Major Minor Major Minor Major Minor Major Minor | 
E W N S E Ww N S E Ww N S E Ww N S| 

0’ Soll! foul’ Saiki Seats 8.29 8.29 8.29 8.29 6.38 6.38 6.38 6.38 4.12: 4.12 eae ane 
5 6.09 8.04 7.61 4.96 Seow fl Oey Hass 3252 4°75)) (3:52) 3.94 2.35) 2°85, Sinton stacy 
10 4.04 7.12 7.00 1.94 - ee AA GION. Pe 223220" 525600 2208 1.69 2.90 1.68 1.0¢ 
15 3.85 5.89 3.94 1.61 DixL Led OSM SRLS 1.48 2.91 1.38 1.35 1.31 260) RES eee 
20 2:39 A507) 2250) 125 OZ Sei Ona teak) 1.02 1.82 1.15 0.80 0.94 1.69 0.78 0.98 
25 > Path. Wek OsO6 1200 2S te 00 00) 0.75 1.09 0.78 0.65 0:69 1.00 “O255) 0248 
30 1 OSs TAG SOR Sill 0.88 1.29 0.72 0.80 0.88 0.95 0.56 0.60 0.56 0.44 0.44 0.60 
35 Ong iesil Os O50 Coe tials Mars O70) 0.69 0.69 0.56 0.45 0.45 0.28 0.31 0.44 
40 0.59 0.75 0.48 0.38 0.55 0.85 0.60 0.60 0.56 0.72 0.44 0.48 0.54 0.45 0.28 
45 OS O88) 0.55 0.75 0.55 0.68 0.69 0.50 0.52 0.45 0.58 
50 0.49 0.75 0.50 0.40 0.55 0.72 0.45 0.60 0.65 
55 0.33 0.38 0.42 0.52 0.60 0.55 0.52 0.58 > 0.48 
60 0.71 0.54 0.48 0.70 0.52 0.78 0.68 
65 0.62 0.40 0.60 
70 0.62 0.46 0.44 


* Densities are given in units of stars per square minute. 
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Fic. 7. Profiles for the major p 
axis. The density p is in units of 
stars per square minute of arc; r is 
minutes of arc. 
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7| Major Axis 


5094 Major Axis 


5101 Major Axis 


5089 Major Axis 


unknown, but probably is of the order of 10” years, 
since its brightest stars resemble those of globular 
clusters. The relaxation time for stellar encounters is 
surely very long, but the smoothing out of large-scale 
density irregularities results from mixing, the charac- 
‘teristic time for which is much shorter. The oscillation 
time for a star through a galaxy of 107 solar masses in a 
700 pc radius is of the order of 10% to 109 years, less than 
the age of the galaxy. 
Another possibility is that there exists in the central 
parts of the galaxy a thin patch of dark material which 
obscures the point of highest density. This suggestion 
“might be considered borne out by the fact that the 
‘deepest plate, ADH 71, shows the asymmetry least 
conspicuously. As has been pointed out, the central 
counts for this plate must be considered with reserve. 
Nevertheless, a dust lane such as that conspicuous near 
the center of NGC 185 (Baade 1944), but weaker and 
smaller, might account for the present observations. 
__ Asimilar asymmetry along the minor axis is derived 
: from the present data (Fig. 10), but is less marked than 
: that for the major axis. 


IV. COMPARISON WITH GIANT ELLIPTICAL GALAXIES 
Mean Profile 


It is of considerable interest to compare the stellar 
density profile of the Fornax dwarf galaxy with the 
luminosity profiles of giant ellipticals. The two can be 
directly compared if the luminosity function of the 
Fornax. dwarf is invariant with radius, a reasonable 
assumption in view of the close agreement shown in 


Figs. 7 and 8. A mean profile for the Fornax dwarf has 
been constructed by smoothing the average major-axis 
profile of plate ADH 5094, according to the other plates’ 
counts and the outer profiles of Fig. 9. The derived 
mean major-axis density distribution is given in Table 
Til. 


Hubble’s Law 


Hubble (1930) found that for 15 normal (giant) 
elliptical galaxies, the distribution of luminosity with 
radius could be represented well by the interpolation 


r 


Fic. 8. Mean profiles for the major axis. Units are as for Fig. 7. 
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_ Fic. 9. Outer profiles from east-west counts. Crosses are E, 
circles W. Lower figure shows schematic W profiles. ( is in stars 
per square minute of arc, 7 is in minutes of arc.) 


formula. 


T=Io/ (r/a+1)? (2) 


logl =logI)—2 log(r/a+1), (3) 


where J is the central intensity and a is a scale param- 
eter. 


or 


Fic. 10. Asymmetry in the Fornax galaxy. The ordinate is the 
ratio of the eastern to the western end of the major axis, and the 
ratio of the southern to the northern end of the minor axis. The 
abscissa is the distance in minutes of arc along the axis. 
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TaBLE III. Major-axis profiles. 


1 5101 5094 5089 : 
a Stars per square minute mean® | 
0’ 7.41 Hots 5.386 3.62 586 
5 6.36 4.91 3.62 2.10 362 
10 4.88 3.66 2.14 1.80 240 
15 4.17 2.89 1.68 1.46 160 
20 2x08 2.10 0.90 0.81 98 
25 Isl 1.10 0.40 0.34 56 
30 0.56 0.52 0.40 0.00 27 
35 0.32 0.34 0.17 —0.14 14 - 
40 —0.03 0.14 0.12 0.00 7 
45 0.12 0.09 0.08 —0.05 1 
50 0.00 =O0-11 —0.07 
55 —0.05 —0.09 0.00 ae 
60 0.03 0.00 in 
65 — Oe? i 
70 —0.06 


8 Member stars per 100 square minutes, based on plate 5094 and smoothed t 
with reference to the others, \ 


The star counts of Fornax can be compared with the | 
analogous formula 


N=No/(r/a+1), (| 


where in the case of this low-density gradient galaxy. 
No can be determined, so that the proper value of a | 
may be derived directly and the formula compared ¢ 


q 
i 
(| 
i 


Computed 


Observed 


log (r+a) 


tS 
Fic. 11. A comparison of the observed distribution of stars — 
(along the major axis) in the Fornax dwarf (crosses and solid line) — 
with that computed from Hubble’s interpolation formula. WV is 
100 times the density per square minute of arc (p). 


with the observations at once. Figure 11 is a plot of 
log N against log(r+a), which for the Hubble law | 
produces a straight line. The value of a is calculated — 
to be 16.2’. As Fig. 11 shows, the agreement with — 
Hubble’s formula is not unreasonable out to about — 
r/a=1, beyond which the observations fall far short 
of the predicted densities. This can be contrasted with 
the case for giant ellipticals, which Hubble found to fit 
the interpolation formula out to values at least r/a= 20. _ 
The fit of Fig. 11 is not improved by choosing other 
values for the parameters No or a, so that we must 
conclude that the star distribution of the Fornax dwarf 
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‘galaxy agrees with that predicted by the Hubble law 


I>) 
yver a much smaller range than in the case of giant 
lliptical galaxies. The reason for this departure will 


be discussed in Sec. V. 


| de Vaucouleurs’ Law 


| 


_ A somewhat different interpolation formula has been 
proposed by de Vaucouleurs (1948), who measured 
three elliptical galaxies and suggested that the following 
formula represented the data well: 


log Z/Io) = —3.33{ (r/a)'— 1}, (S) 


where.a is the radius within which half the light of the 
galaxy is emitted, and Jo is the luminosity at a. He has 
since suggested that all ellipticals so far measured fit 
this formula better than Hubble’s (de Vaucouleurs 
1959). 

_ To test the Fornax galaxy’s agreement with this law, 
we have in Fig. 12 a plot of logN against r*, which for 


log N 


Fic. 12. The distribution of stars, plotted in units in which de 
Vaucouleur’s luminosity distribution law makes a straight line. 


de Vaucouleur’s formula gives a straight line. The figure 
shows that there is no agreement whatever with the 
formula, within the distances to which observations can 
be carried. 


V. GRAVITATIONALLY LIMITED RADIUS 


As shown by von Hoerner (1957), there is a limiting 
radius for star clusters in the vicinity of a large mass 
(the Galaxy) beyond which members are lost to the 
cluster. It is at the point where the acceleration on the 
star due to the Galaxy is balanced by that due to the 
cluster. The expression for the gravitationally limited 
radius is (approximately) 


r,=R,(m/2M,)*, ; (6) 


where R, is the distance from the-star to the Galaxy, 
M,, is the mass of the Galaxy, and m is the mass of the 
cluster. The radius thus computed has been compared 
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Fic. 13. The distribution of stars plotted against the reciprocal of 
the radius, to determine the limiting radius. 


with the observed radius of M15 by King (1961), and 
with the sizes of intergalactic tramp globular clusters by 
Burbidge and Sandage (1958), who point out that the 
appropriate value for R, is the distance of closest 
approach of the object to the Galaxy, rather than 
necessarily the present distance. 

As Fig. 11 demonstrates, the star density in the 
Fornax galaxy falls off rather abruptly compared with 
the Hubble law, and it is tempting to believe that this 
is the result of the gravitational influence of the Galaxy 
on the outer stars. We may use Eq. (6), assuming 
reasonable values for the mass of the galaxy and its 
distance, to derive an approximate value for the gravi- 
tationally limited radius of the Fornax dwarf to compare 
with the observed radius. Similarly, we may turn the 
problem around to see if its observed radius confirms 
our rough estimate of its mass. 

Adopting 188 kpc as a reasonable value for the 
distance of the Fornax galaxy from the Galactic center, 
and 1.8X10" Mo as a reasonable value for the mass 
of the Galaxy (Brandt 1960), we have left to estimate 
the mass of the Fornax dwarf. 

One method is to compare the observed portion of 
the luminosity function of the galaxy, as determined 
from the star counts, with that for M3, for which 
Sandage has estimated the total mass (1954). The total 


ite) 
No. 
5 
g R 50’ 100° 
Distance 


Fic. 14. The distribution of background galaxies 
behind the Fornax dwarf. 
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Taste IV. Galaxies in the field of the Fornax dwarf. 
Number of 
Number of galaxies 
Radius Area of ring galaxies per 1000 
of ring (square minutes) in ring sq minutes 
10’ 314 2 6 
20’ 943 5 5 
30’ 1570 13 8 
40’ 2200 11 5 
50’ 2827 12 4 
60’ 3456 18 5 
70’ 1000 (sector) 6 6 
80’ 1000 (sector) 4 4 
90’ 1000 (sector) 4 4 
100’ 1000 (sector) 4 4 


number of stars (which are members of the galaxy) on 
plate 5101 is found to be 4295. The difference in mag- 
nitude between the brightest stars in the galaxy and 
the limit of the counts is not known accurately, but is 
estimated from measures on photoelectric sequences on 
other ADH plates exposed the same night and for the 
same length of time to be approximately 0.7 mag. (pg). 
For M3, Sandage (1954) finds ~50 stars in the top 
0.7 mag. Thus, the comparison suggests that the galaxy 
has approximately 85 times the number of stars in this 
interval as does the globular cluster M3. If we assume 
the luminosity functions are the same, then we may 
derive a mass for Fornax by multiplying the mass of 
M3 (Sandage 1954) by this factor, giving 


m=1.2X10' Mo. 


It should be pointed out that, although the bright stars 
of Fornax do resemble globular cluster stars, there is 
the suggestion that the luminosity function of the 
galaxy is not identical to that of a typical globular 
cluster, either at the bright end (Hodge 1961) or the 
faint end (Sandage 1954). 

Another estimate of the mass of Fornax, also rather 
uncertain, can be made from its total apparent mag- 
nitude, assuming a mass-to-light ratio. Unfortunately, 
both of these quantities are poorly known. The total 
photographic magnitude has been estimated roughly by 
eye to be approximately 9.0 (Shapley 1939). The mass- 
to-light ratios for giant ellipticals are found by Page 
(1959) to average approximately 30, using a Hubble 
constant of 75 km/sec/Mpc (Sandage 1958). However, 
it is probable that dwarf galaxies more nearly resemble 
globular clusters, which have much smaller ratios. If 
we adopt a mass-to-light ratio of 2 (Sandage 1954), the 
mass of Fornax is estimated to be 


m=2X10' Mo, 


in surprisingly good agreement with that derived from 
the comparison with M3. 

Proceeding with an assumed mass of 2X10? Mo, we 
derive a value of the tidal limiting radius, r,, of 7.0 kpc. 
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The observed radius may be determined using the} 
method employed by King (1961). Figure 13 shows al 
plot of the reciprocal of the radius against V? (where Vj 
is the projected density of stars in units of number per} 
square parsec). The radius obtained in this way is 56’,} 
which is 3.3 kpc, assuming the distance of 188 kpc. Thus} 
the observed limiting radius is approximately half the} 
computed tidal limit, a very reasonable result, sug-} 
gesting that at some time Fornax has approached the} 
Galaxy somewhat nearer than it now lies. Its present] 
radial velocity relative to the center of the Galaxy is! 
—74 km/sec (Humason, Mayall, and Sandage 1956), } 
which amounts to 74 kpc per 10° years, so that a closer? 
approach during history is a very reasonable possibility. } 
i 

VI. TRANSPARENCY OF THE SYSTEM al 

In order to check the transparency of the Fornax } 
dwarf galaxy, to see if general internal absorption | 
might affect the counts, a survey of background galaxies | 
seen through the system has been made. To be sure of } 
identification, only galaxies with decidedly elongated | 
images were noted. All such objects were counted in | 
circular rings of 10’ width, centered on the Fornax | 
system. The results are given in Table V and shown in 
Fig. 14. As these show, the frequency of galaxies behind | 
the Fornax dwarf (radius 56’) does not differ from that | 
in the surrounding sky; there is probably no significant | 
amount of general absorbing material in the galaxy. _ 


VII. SUMMARY 


The study of the distribution of the brightest stars - 
of the Fornax dwarf galaxy has led to the following - 
conclusions: (1) The gradient of stellar density (and 
luminosity density) is, exceedingly small, compared to 
that for giant elliptical galaxies. (2) The stars in the 
Fornax dwarf are distributed quite smoothly, with no- 
real irregularities detectable. (3) The major axis is 
oriented so that it lies 54° eastward from the north. 
(4) The ellipticity of the equidensity contours varies 
with distance, increasing from 0.21 near the center to | 
0.36 in the outer portions. (5) The distribution of stars — 
near the center of the galaxy is asymmetrical; possibly 
a small, thin patch of absorbing material obscures part 
of the center. (6) The major-axis profile agrees reason- 
ably well with the Hubble law out to r/a=1, beyond 
which the galaxy’s stellar distribution falls off much 
more rapidly than the Hubble law. (7) There is no 
apparent agreement with de Vaucouleur’s interpolation 
formula. (8) The abrupt cutoff of stars in the galaxy 
can be explained on the basis of the tidal limit on its 
radius imposed by the mass of the Galaxy. (9) The 
observed limiting radius of the Fornax dwarf is 3.3 kpc, 
assuming Baade and Hubble’s (1939) distance of 188 _ 
kpc. (10) A check on the distribution of background 
galaxies indicates that there is no large-scale general 
absorption within the Fornax galaxy. 
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The Motion of a Hyperbolic Artificial Satellite around the Oblate Earth 


GEN-ICHIRO Horr 
Yale University Observatory, New Haven, Connecticut 


(Received May 15, 1961) 


The first-order perturbation of the oblateness of the earth on the hyperbolic motion of an artificial satellite 
is developed by the von Zeipel method. The theory is valid for any eccentricity greater than unity and for 
any inclination. The modification of the Delaunay variables and the meaning of integration constants are 


discussed in some detail. 


1. DELAUNAY VARIABLES IN HYPERBOLIC MOTION 


F a stands for the distance of the perigee from the 
center of the hyperbola, and ae the distance between 
the focus and the center, 


r=a(e—1)/(1+e cosf), (1) 


where f is the true anomaly. Furthermore, 7 is defined 
by Kepler’s third law, w=n?a’, and / by J=nt-+ const. 
For convenience we call a the semimajor axis, m the 
mean motion, and / the mean anomaly. The notation 
of J for the inclination, # for the longitude of the 
ascending node, g for the argument of the perigee, and 
e for the eccentricity, corresponds to that in elliptic 
motion. 

Let c;(j=1, ---, 6) be a set of osculating elements. 
The equations of motion in these variables are 


6 OR 
= [ex,6; Jé;= ? k= i‘ ass 6 
j=l 0 


Ck 


where R is the disturbing function, and [c;, c;] the 
Lagrange brackets, 


A(x,%) A(y,y) (2,4) 
a(p,q) 9(p,9) 9(p,9) 


P, q being any pair of cx, c;. If (@’, y’), (@”, y”), and 
(x"’’, y’”) are new coordinate systems obtained suc- 
cessively by rotations by the angle h about the.z axis, 
the angle J about the x’ axis, and, finally, the angle g 
about the z’’ axis, we have 


[¢,q]= 


a(h, PN eye) 
Pere eee 
3(,9) 
CZ, Ti3. (hae VE FP 
(bq V Eee ee 
8(p,9q) 
s a O(g, x Jt} yl —a"yl") 
Leal’ =Loq"+ 
(p,q) 
or 
G) 8(h,H 
Celso ee 
(p,q) a(p,9) 


if G and H are defined as in the elliptic case by 
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G= ay aly AA ies = =([pa(e— 
H=Gcosl. 


1)}, 


x” = a(e—1), y'"=0, 
#/”"=0, in Cts 
y= ¢ —) : 
oe m ae—1 
a?(e—1)? =o 
Then we find 
Oca) 
qi'=>) 
3(p,9) 
poo O(e,k) O(a!” ,4'"") 0(k,a) 
d(e,k) (p,q)  A(k,a) O(p,g) 
; O(a!" a") aa,2) | 3 
a(aje) d(p,9)4 
atk) (a,e) 9(p,9) | 
0(p,9) 
if 
l=n(t+k), 
and 
C=) 2a: 


The canonical equations for C, k are 
C=0R/dk, k=—dR/0C. 


If we use / in place of k, 


OR 10R | Gu ( 
al. yy n 2 


Van “f (ua)"] 
dt ¢ 


a 


which leads to the definition of L: 


L=-—(ua)}. 
Finally, 
5 OK Ol ane oF 
j=n-—=-_(—-R) =—— 
dE JLN2E oL 
if ‘ 
— (i /2L)4R. | 


In summary, the relations between the Delaunay 


/variables and Kepler’s elements are 


=— (ua), 1=nt+const, 
G=[wa(e?—1)}?, g=argument of the perigee, (2) 
H=G cosl, h=longitude of the node. 


] 
|The equations of motion are 


| dL/dt=dF/al, dl/dt= —8F/aL, 


dG/dt=0F/dg, dg/dt=—0F/dG, (3) 
dH/di=0dF/dh, dh/dt=—0F/0H, 
with the Hamiltonian, 
P= —(ui/2L*)+R. (4) 


In this derivation the method used by Brouwer and 
Clemence (1961) for the case of elliptic motion was 
followed. It was essential to go through the derivation 
in detail in order to ascertain the form of the equations 
in the Delaunay variables for hyperbolic motion. 


2. THE DISTURBING FUNCTION 
The disturbing function is 
R= (uk2/r*) (1—3 sin’B)+O(k:”), 


with the equatorial plane of the earth as the reference 

plane and @ the latitude; n=k?M if M is the mass of 
the earth, & the Gaussian constant. With the Delaunay 
variables, the Hamiltonian has the form 


F=F)+F,+0(k.2), 
| ui Re Ib ONG 
alee 
15 Di ERECT AG 
| SP) UENO 

Hee) Sclictn| 
{ 2 24/r 


which is of the same form as in the elliptic case except 
for the sign of the Keplerian term; a*/r? and (a3/r*) 
Xcos(2g+2f) are understood as functions of the 
Delaunay variables through the Keplerian relations: 


Fo=—p/2L?, 


r=a(ecoshu—1), /=esinhu—u, 


u being a parameter. 


3. CANONICAL TRANSFORMATION 
A canonical transformation, 
(L,G,H,l,g,h) — (L’',G’,H,I’,g',h’), 
is performed with the use of the determining function 
S=L'1+G'g+Hh+S,(L',G',H,l,g,—)+0(k2) 


so as to eliminate an angular variable or variables from 
the Hamiltonian. However, in contrast with the elliptic 
case where a*/r* or (a*/r*) cos(2g+2f) consists of 
periodic and secular parts, there is no such distinction 
in the hyperbolic case. Therefore the choice of the 
determining function and the associated new Hamil- 
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tonian F* is not unique as it is in the elliptic case. 
Different choices may result in different meaning of 
integration constants introduced. The choice must 
depend on simplicity of the development of the theory 
and convenience of application. 

Provisionally we assume 


The energy integral 
a 
will then give 
Fut= Fo(L')=—1t/2L", 
OF OS, 
—+tF,(L,G,A,Lg, re. ) =0. 
OL ol 
The equation for Sj is 
OS, wWkf @& a 
= E +B: cos(2e-+2,) (6) 
al IEEE re Te 


where 
A=-3+4('/6"), 


The new equations of motion are 


B=3-}('/G"). 


L/=0F*/al, U=—oaF*/aL', 
G'=0F*/dg', g'=—dF*/dG’, 
H=0F*/dh', h’=—0F*/0H, 


with the Hamiltonian 
Py oe — ee 


Integration leads to the solution 


= const ant —=(— ES) t-- 105 
G’=const, g’=const, 
H=const, h’=const, 


with L’, G’, H, ly’, g’, h’ as six constants of integration. 


4. INTEGRATION OF 0S;/0l 


Equation (6) is integrated with use of the area 
integral in the form 


a’ Lora Ike 
[ranrne fram = f atecosnar 
im GJ r G 


153 : 
Poa sinf), 


a 
— cos(2g+2f)dl 
(es graf 


1K: 
tae fate cosf) cos(2g+2f)df 


L 


é é 
--|- sin(2g+2f)+— sin(2g+f)+—- sin(2g+3) | 
G*L2 2 6 
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and finally we have 


uke it 
SiO | 4(ftesing) +B] sin Qe+2h 


oe sin(2g+f ES sin(2g+3f ] | 


+2Z(L',G',H,g), (7) 


where Z is an arbitrary function, the introduction of 
which is not needed in the elliptic case in which Sj is 
determined by the periodic part of F;. The determina- 
tion of Z is necessarily associated with the meaning of 
the constants of integration L’, G’, H, ly’, g’, h’. 

The relation between Z and L’ is 


L=L'+(aS;/él) 


2 a? 
= hf = E —+B— =, cos(2e-+21) | 


which shows that L’ is the value of L at r= or 
t=-+«, The meaning of the other constants should 
be consistent with that of L’. In order to make the 
problem completely defined, a choice must be made 


between imposing the conditions for = — © orf=+ 0. 
If the former is chosen, we shall have, for ‘= — ~, the 
conditions 

G'=G, =I, g =, h'=h, (8) 


which will determine Z uniquely. 


5. DETERMINATION OF Z(L’,G’,H,g) 
Equation (7) leads to the relation between G and G’: 
G=G'+ (dS1/dg) 


2 


iT 
= 6 +9 cos( 25-49 )+e cos(2g+f) 
aZ 
+- cos(2e+3f) [+ 
3 ae 


Assuming the direct motion of the satellite, we have, 
at f[=— , 


1 (e—1)? 
cosf=—-, sinf=— , 
€ € 
De 2(e—1)? 
cos2f= sin2f=—_—__— 
e? e 
3e—4 (2&—4)(e&—1)? 
cos3f= , sin3f= ; 
a 


and 


cos(2g+ 2f) +e cos(2g+f) - cos(2g+3f) 


GEN-ICHIRO HORI 


From Eq. (8) it now follows that 


Zs whe [3-2 7(@—1)m 
=— a Ye sn2g| 
e e 


and integration leads to 


sin2g— 


3—2¢ 3(e—1)3 
1 cos | 
G’ 


} 
i 
©) 
where Z; is again an arbitrary function. Substitution oi | 
Eq. (9) in Eq. (7)_results in 


4-2 (LGB), 


2B» 


=| | 4(tesinn) +3] - sin(2g+2f) 


Si= 


e e 
te sine sin(2g+3f) 


He-1)h 
5 ae cost | 
e 


+2Z,(L’,G’,H). 


(7’)| 
6. DETERMINATION OF Z;(L’,G’,H) 


The arbitrary function Z, will be determined by the | 
conditions (8) after expressions for J, g, and h are’ 
obtained. Differentiation of S:, Eq. (7), with respect — 
to L’, G’ is performed with the aid of the relations: __ 


ov (e) 1@oy dle) 1G dy | 
oL a eLide aG oe L de | 

0 fa a of On AE? | 
— -) ==> Os, == = (“+-) sinf. | 
de\r fe de r @ mi 


The result is for /, 


V=P (6S 00) 


1 wk @G? «a 
=/+— {4 ree Lo 
e L2G’ PL? ¢ 
+3]( a@G? a ) ( 
—B) {| —— —~—-— 1 } sin(2g+ 
2 Pl ety 
—— ) 4 
— — -—-— } sin(2¢+3f) —— sin2: 
PrP Lr 3 Ses 3e 
4+ 2 G’ OZ, 
— costg]| | +—. 
, OL’ 


The coefficient of $B in the right-hand member should 
vanish at t= — ~, since the corresponding coefficient of 


ee 
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Bin Eq. (7’) is identically zero as a function of LOG 
H, g at i=—. Z; is, therefore, determined by 


Var whe (@=-1) wkoA 
Sere 6 8 1(L2+G") 
_ After integration, we have 
jin fi 1 i 
a a i —+— tant )+2.(6'H) 
G? eG! G! 


wke 1 

= 4 @—pH+tan ——] 4.2.16, (10) 
-G’ (e2—1)3 

_or substitution in Eq. (7’) gives, for Sj, 


‘ wR 
Sy= 
} G8 


1 
{4 [ie sinf+ (e— Detar] 
e—1)t 


1 
+5|- sin(2eb9f)+- sinQe+f)+~ Bay +3y) 


1 ek 
+(--— 
2 Ses 


+2Z2(G',H). 


oa 1)3 cos | 
3e 
(7"’) 


7. DETERMINATION OF Z»(G’,H) 


Differentiation of Eq. (7’) with respect to G’ gives 
_ the expression for g: 
—| 
€ 


- g’=g+(dS1/0G’) 


ECS 
a ae 1, | sinf— 
TG? ele + 


1 2G? a 
‘en 
2 PL? ¢ 


a@G? a i 4 
+(—t-—-) sine +3))—— sin2g 


PL” 


A+-2¢ G! hs 
dlr ey 
3 OL! Gis 


1 
x| He sinf+(e—1)}+tan va | 


@2@G? a 


1) sin(2g+) 


Si Lou * 
2 a G? 


(- -—)|; in(2g+2f) + sin(2g-+) 
me os sin(2g-+ 2f Boner ts 


1 
+ 2e+3 +(-.) 2 
| sin(2g+3f) aa sin2g 


ee 1)3 cost | 
3e 


The coefficient of B in Eq. (7”) is identically zero at 


=—o, hence the coefficients of 3B and 
9/2—(15/2)H?/G”? 
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in the above expression vanish at t=— ©. The coef- 
ficient of A in Eq. (7”) is, on the other hand, expected 
to be identically zero or at most constant as a function 
of L’, G’ at t=— ~~. In fact, we see that 


1 
Ire sinf+ (e—1)?+tan 
= 1)? t=—00 


Therefore the coefficients of A and —3$+(15/2)H?/G” 
in the above expression are 0 and $7, respectively. 
The equation for Z» is then 


3 whe(_3 15 HP 
Ge) 
o GaN a2 


which by integration yields 


022 
aG’ 


Zo= 


3) Hay /ll Sse 
or (-—-—-) +28) 
DEGEN G2 


2 


3 uke 
= a Aa) 
2G 


Substituting in Eq. (7”), 


wR 1 3 
S1= {al H+< sinf+ (e—1)?+tan ->] 
G8 (e—1)? 2 


1 
+5] sin(2g+2))+- sin(2e-t/) 


1 
* sin(2g+3f)+ (---) sin2g 
a 


(eeiBy 
a 


cos | }+Z0(H) 


3e 


The coefficients of A and B are now identically zero 
at t=— and Zs;, an arbitrary function, is at most a 
constant which can be put equal to zero. The final 
espression for S; thus reduces to 


Whe if saya: 
Si= | (-345 —) fe sinf+(e?—1)} 
G Dea? 
3 Seo Ae Nip 
+tan“ -=|+(-->=) sin(2g+2f) 
(Ak! 2 2\52'G2 


1 
i sin(2g+f)+— * sin@agt3)+(———) sin2g 
e 


1 
+-——(e—1)! cos | |. (11) 
Je? 
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8. SUMMARY OF THE FIRST-ORDER PERTURBATIONS 


whos 1 3H \a 
(Le) 
L’ 2 2G?/r 
3 Sse 
+(---=)= cos(2g’ +2) 
2 2G"? 
whe 3 3 
GG ce (=) costoe' tap 
CONDO GZ 


€ 
+e cos(2g’+f) ne cos(2g’+3f) 


2 2G? 
+(1-=) cos2g’ oe — sin2dg ‘| 
Se gids 
H=H, 
1 whe (7 1 3 
jnr--—*|(—-+-—) 
CLE 2 @? 
a G 7a (e—1)3 
x|(5 Wie i) sinf— | 
ye Le e 
Sak OG? “@ 
+(--- YI —------— 1) sin(2 +f) 
4 4G? je Ie 


AN Oa 


+ = 
yr L” 


Cou 
+==) sin(2g’+3f) 
tas 


4 4+ 2 G’ 
—— sin2g’+ — cos2g 
3e Se ae 


1 wk» 1 
es Sorc {( ce 
Ch AG" 2 


| il 


2 a, 


2G? a G’ 
Spel ) sinf-+—| 
pL? r eL’ 


3 3H eG? o 
4 4G" Pie ¢ 


x 


i) sin(22’+f/) 
eC as 
yr Sa 


=f ) sin(29’+3/) 


4 4+ 2 G’ 
——— sin2e' = — cost | 
3e Seeks 
WR ou As =) G’ 
a (+ ) (He sinf—— 
G’ yy 2 GE ib 


LE 3 )+(- - 1 ] 
—tan ——-_7 ——— — }]-sin(2e’+2 
G 2 2 QeG? ; one, 


e e 
- sin Cen ome sin (2¢’+3f) 
eal "3 
a 
2 3é 


sin2g’ -—— a cos2g 


|} 
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(Ge 
--é sinj——-— taney aaa 
ibe Gee 


Suk 


aa 
h=h— oA als 


1 Cat e 
ie sin (2g’+2f) ns sin(2g’+f)— ‘ sin (2g’+3f) 


abe sel 
_ (--—) sin2g’ 
De waves 


{Gta 


1G 
+— — cos2g’ 
Semi 


where 


and a/r, f are the functions of e, /’ through the implicit” 


relation with u as a parameter, 


a 1 e—coshu 


>] 
r ecoshu—1 


) 
e coshu—1 
l’=e sinhu—u. 
The new mean anomaly /’ is a linear function of the time 


fa(— pope, 


and Jy’, together with L’, G’, H, g’, h’, constitute six 


integration constants. 


9. CONCLUDING REMARKS 


The development of the theory is quite analogous to 


(12) 


the elliptic case treated by Brouwer (1959). The essen- 


tial difference is, however, seen in the meaning of the 


integration constants associated with the choice of the — 


new Hamiltonian, Eq. (5). With the use of the von 


Zeipel method, this choice may be the most suitable, 


and probably unique one, for the hyperbolic case. 


The essential feature of the present first-order theory — 


may hold if the solution is extended to the second and 


higher orders. We may adopt the same choice for the 
new Hamiltonian and then the transformed Delaunay — 


variables are constants or the same linear functions of ~ 


the time, with the same meaning as in the first-order 


theory. In addition, as there is no separation into — 


short-period, long-period, and secular terms, we have 
no critical inclination. As we cannot expect, therefore, 


any new theoretical feature in the higher-order theory, — 


and as the perturbations decrease after the satellite 
passes its perigee, the present first-order theory may be 
enough both for theoretical and applicational view 
points. 


10. MODIFICATIONS REQUIRED IN THE 
REPULSIVE CASE 


The present theory was first developed for the case 
of repulsive force (u<0) in an application to the scat- 
tering problem. The results are quite similar to the 
present case of attractive force. The following is a 
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6 —|( ale —)[ ; 
= 7 Ga Slil 
edge 


summary of formulas that apply to the repulsive case: 


r=a(e—1)/(e cosf—1), 


l=e sinhu+u, 


r=a(e coshu-+1), 
—p=na'; 
Delaunay variables: 
a= — [(—p)a}}, 
G=[(—n)a(e—1)}}, g=argument of the pericenter, 


H=((—n)a(e—1) }* cosf, 
h=longitude of the node; 


l=nt+ const, 


Hamiltonian: P=— (u?/2L?)+R; 


Disturbing function: 


uke IL SFE NCE 
LPN DO Coe 


Si Orble Wa? 
+(---=)= cos(2e+2p) | 
2 2@/r 


Determining function: 
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1 7 3) Sek 
—(?—1)}— tan“ = }+( — ) 
(2 een 


1 e e 
x| 5 SUNG St ia a2 al sin(2g+3f) 


HE 1 ) - (e—1)3 
S= == |) SIMe— cos2 || 
Die: 3e? ° 


Comparison with Eq. (11) for the attractive case shows 
various changes in sign, and the change from —3z to 
+4 in the coefficient of the first part of Si. 
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The authors present an approximation to exp(—6) in powers of 6 that is more convergent than a Taylor’s 


series development if @ is limited to a given range. 


N the development of our theory (Brouwer and 

Hori 1961), we used the Taylor’s expansion for 
exp(—aax), where a is the coefficient of the radius 
vector in the exponential model of the atmosphere, a 
the semimajor axis, and *x=(r/a)—1. Since the ex- 
pansion is not effective when |aaa| takes values greater 
than unity, we suggested in the paper to use the approxi- 
mation of the form 


where the coefficients A; may be determined by the 
principle of the least square when a range of @ is given, 
Say, —yS0 < aims 
For simplicity let C; be 
Ci= (—1)'(4;5/9)) 


the evaluation of C; is then formulated by the equation 


7=9, 1, “oe 


Y 5 
i= f ey e%—S* C,67?dd=minimum, (1) 
j=0 


= 


where expL— (6+) _] is a provisional weight function 
taking the value 1 at the perigee (@= —7) and exp(—2y) 
at the apogee (@=+7). 
The conditions 
61 /eC=0, 


lead to the normal equations for C; 


g=0 Les 


b] 


y 


fe 
Qitke- C+) G9= f Gke— C81) dp, 
ea ea 

B=0,1,-+*,5: (2) 


Evaluation of the integrals results in 


5 
DG; 
j=0 
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it wen dom | er = (eee | 

se 274 k=0 R! k=0 R! 
7=9, ue fay 


and 


if 


head! a 
viedo] —(—y)'-e27 > I 
k=0 k! k=0 k! 


j=0,1,--- 


li Q;, Q;' (j=0, 1, ;-++, 5), and Pj, P; (G7= Onis. 
are defined by | 


Qo=l, QF=QiHrt (— 127/79), : 
Qo'=1, Of =O;-1'+ (2/797; 7=1,2, --755, | Ge 
Po=1,  Py=Pyat(—D(1/)y’, | 
PJ =1,. P/ =P + 0/73 7=1, 2) ee 


and b;, Qjk by 


a 
-— —_(() .ey—(),'¢-3 _ wee 
b= (et Oe), J 0,1, +5, 


a;= j\(P;—P,'e”); j=0, 1, RO (4) 


j, k=0, 1; os 


Qjk= Qk, 
the normal equations are 
5 


DX ajxCn= by, 


k=0 


j=0, Ue IES, (5): 


The formulas (3), (4) show a method of evaluation of 
the integrals appearing in the normal equations (2). 
The normal equations (5) were solved for five cases, 
y=1,2, ---,5 with the use of the IBM 650 and IBM 
610. The resulting values of A; are given in Table I. 


TABLE I. The coefficients A’; for y=1, 2, ---, 5 and their mean errors. 

y 1 2 3 4 5 
Ao +0.9999900 +0 .9972809 +0.9927548 +0.7834172 —1.0892184 
A; +0 .9966676 +1.0143125 +1.0874261 +1.5197473 +2.8373974 
A» +1.0002082 +1.0054912 +0 .9636549 +1. 1811823 +3. 2873689 
A; +1.0840814 +0.8701602 +0. 6360486 —0.4107438 —3.4678006 
Ag +1.0209843 +1.0931240 +1.3336037 +1.3402296 +0. 2884029 
A; —0.3977536 +1.7402653 +2.0209304 +3 .2974657 +6.0651410 
in 0.00058 0.0039 0.018 0.12 0.55 

Ho 0.0006 0.027 0.35 253 11 
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DONG SOUUTION FOR 


The mean error p is evaluated by 
{ f= /(i—<"), 


where J is the integral defined in Eq. (1). 
J If we put 


Boy =C7?+2 >D. (Cosh Crane j= 0, ile vee, 5, 
k=1 
Boj41=2 »y One Garaena j=), ile 29D 4, 
k=0 


the C’s under the summation being limited to Co, Ci, 
+++, Cs, we have 


ne 


| oY 5 
/ if e~ 84+) da— 2 s C; O5e— 28+7) de 
iF j=0 Jy 
10 a 
+> B; | o%e-¢+ a0 


0 J_y 


5 10 
=1(ey—e-47)—2 D Cyb; + Bja;. 
j=0 ee 
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For a comparison, the mean error po is calculated in 
the case of the Taylor’s expansion (A4;=1, j=0, 
il * cin 5) : 


1 Y 
Loe=— 
dy J_y 


iene Niger 
|e (1-04 P—-4+—H— )| d6. 
Dy oe 24130 


The results for u, “o are shown in Table I for the five 
Casesay. =e Dien 5s 


ACKNOWLEDGMENTS 


This investigation was supported by the National 
Aeronautics and Space Administration. It is a pleasure 
to acknowledge our indebtedness to Dr. Y. Kozai. 
Discussions with him on the possible divergence of 
developments that arise in the drag theory of artificial 
satellites led to the procedure developed in this paper. 
We also extend our thanks to Mr. W. Wilkinson, staff 
member of the Yale University Computing Center, for 
aid with programming on the IBM 650. 


THE ASTRONOMICAL JOURNAL 


VOLUME 66, NUMBER 6 CIGU Sie 


The Orbit of Comet Schwassmann-Wachmann 1 
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A numerical integration and orbit correction for Comet Schwassmann-Wachmann 1925 IT reveals the 
existence of nongravitational forces acting upon the comet. The comet approaches within 1 a.u. of Jupiter 
in 1974, and drastic changes in the elements will be produced. Extended opposition ephemerides are 


furnished until 1973. 


HE Comet Schwassmann-Wachmann 1925 II is a 

uniquely peculiar object in the solar system 
because of the sporadic large increases in brightness 
which occur without any known explanation. To provide 
the observers with adequate ephemerides, a set of 
orbital elements was established (Herget 1947), based 
on observations extending from 1934 to 1944. Sub- 
sequent ephemerides have been published on Harvard 
Announcement Cards. An extension of this work has 
now been undertaken by means of the standard IBM 
650 electronic computer programs which are in regular 
use at the University of Cincinnati for the determination 
of minor planet orbits. : 

Nearly 200 observations of the comet were selected 
haphazardly from the literature with the view to cover- 
ing rather uniformly the entire observed arc. The 
resulting orbit correction gave residuals in right ascen- 
sion which are shown in the upper part of Fig. 1. There 
are systematic deviations in excess of ten seconds of 
arc. Every plotted point was included in the least- 
squares solution, although several outstanding points 
are of doubtful validity. The declination residuals are 
not shown. 

A second solution was then computed, excluding 
those points which are marked off by a dashed curve, 
as shown in the lower part of Fig. 1. The residuals of 
this solution are shown in Table I. The amplitude of 
the systematic deviations is reduced to about three 


+10~ 3 4 


seconds of arc, and there appears to be a sharp discon- 
tinuity near JD 2426000. It may be assumed for the 
present that lesser discontinuities, which have not been! 
resolved, exist at. other points. If we adopt a residual 
of 20” in the mean anomaly at JD 2425200, this is’ 
equivalent to a change of 07025 in the daily mean 
motion, corresponding to an impulsive change in the 
velocity of the comet of one meter per second. The 
previous orbit (Herget 1947) had not shown any evi- 
dence of such deviation from gravitational motion, but 
it may be seen from Fig. 1 that this was entirely 
fortuitous. | 

In the future, the comet approaches Jupiter within 
2.0 a.u. during the interval from JD 2440960=1971 
Jan. 9 to JD 2443600= 1978 April 2, with a minimum 
of less than 1.0 a.u. near JD 2442220=1974 June 22. 
The nature of the changes in the osculating elements 
may be seen from Table II. After the close approach to 
Jupiter, the eccentricity is diminished to about half its 
present value, and the period is diminished to such an 
extent that there is not another close approach to 
Jupiter for nearly a century. It is not intended to imply 
that these are accurate prediction elements, since other 
sporadic changes may occur, and the perturbation 
computations are not entirely adequate at the 20-day 
interval which has been used. However, the printed 
ephemerides (Table III) may be expected to be reason- 
ably reliable for practical purposes. 
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aE cosd Ae 
5 
5 
5 
8 
7 
9 
6 
5 
9 
0 
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6 
7 
6 
7 
6 
7 
9 
9 
9 
8 
7 
8 
8 
8 
5 
7 
7 
7 
9 
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9 
0 
7 
8 
7 
7 
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9 
9 
6 
6 
7 
9 
0 
7 
7 
7 
7 
9 
9 
9 
8 
8 
6 
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9 
8 
7 
6 
7 
0 
9 
0 
8 
8 
9 
8 
9 
7 
7 
7 
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cosé Aa 
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TABLE I (continued) 


JDE—2400000 cos6 Aa = AS cosé Aa Aé Place | JDE—2400000 cosi Aa <Aé6 cosé Aa Ab 
356 66.7 — 4"7 + 078 — 0%9 + 178 Lick 364 31.9" + (001) 027 + 1%4 + 2"1 
356 66.8 —5.4 + 0.5 — 1.7 + 1.5 Lick 3643858 = 1300 2.9 + 4.9 + 4.3 
356 88.7 — 3.7 + 0.4 0.0 + 1.4 Lick 364 58.6 + 0.1 = 2.9 + 1.3 — 1.4 
356 88.8 — 4.9 + 0.8 — 1.1 + 1.7 Lick 364 68.8 — 1.2 + 3.2 +7051 + 4.7 
357 19.6 —69 — 2.4 — 3.4 —) tO") Werk soe Wins Sele ae 28) + 2.3 + 4.3 
360 19.8 +0.6 + 2.3 se ail + 3.4 Yerk 365,19.6) > —"3.2> =-"0.2 — 2.1 + 1.6 
360 43.8 — 4.6 + 6.1 — 2.9 + 7.3 Yerk 365 36.5 + 1.0 =F 16 am ee) + 3.0 
360) 71.8) %=9 353." 2082 — 1.5 +- 1.0% Lick $65, 59'5) 1 OE On — 0.8 189 
360 71.8 —3.8 -—0.1 — 2.1 + 1.1 Lick 365 43.5 — M2 & 2.4 — 0.2 — 1.0 
360) 73.7 | — 31> =S/055 — 1.3 + 1.7 McD 365 46.5 °F 1.8 — 1.3 + 2.8 + 0.1 
SOL 05.7, 2 — 3.2) =e — 1.6 + 2.4 Lick 365 49.5 —0.3 +4 6.0 + 0.7 + 7.3 
361 05.8 — 4.1 +1.7 — 2.4 + 2.8 Lick 36578. - — O15" = 2:3 + 0.4 — 1.1 
361 06.5 — 2.6 + 0.4 — 1.0 + 1.5 Yerk 367 59.8 —5.1 — 3.4 — 3.2 — 1.8 
361 39.6 —3.1 + 1.1 — 1.7 + 2.1 Yerk 368 12.8 — 4.5 — 3.1 — 2.2 — 1.1 
361) 53.5) = 34-5) = e2e — 2.0 — 1.1 Yerk 368 13.8 —16.0 —0.1 (—13.7) + 1.9 
364 27.8 — 2.4 40.7 — 1.1 + 2.1 McD 368 14.7 -—17.0 — 1.9 (—14.6) + 0.1 
364 29.8 +04 +4 2.9 + 1.7 + 4.3 Yerk 368 19.9 —1.4 — 0.6 + 1.0 + 1.5 
364 29.8 + 2.2 + 3.4 + 3.4 + 4.8 Yerk 369 01.6 + 1.8 +11.9 + 4.0 (+14. 1) 
364 30.9 —0.8 +04 40.5 + 1.9 McD $69 32-5. ) 7365 9252 + 9.6 — 0.2 

TABLE II. Orbital elements. 
Epoch JD 2430000.5=1941 Jan 6.0 ET 
Mo 353°57207 (1950.0) Je Q 
n 0.06103951 w 356°22229 +0. 74633139 +0.65774613 
a 6.3884541 Q 322 .00326 —0.60620675 +0 .60863191 
e 0.13550168 4 9.51656 —0.27477780 +0.44377609 
Epoch JD 2436000.5 Epoch JD 2442000.5 
My 1°79256 (1950.0) Mo. 354°89931 (1950.0) 
n 0.06120884 w 355°81213 n 0.06498197 w 13°38059 
a 6.376666 Q 321.60963 a 6.127377 2 319°66031 
€ 0.13156612 1 9 48662 e 0.11077967 1 9.75739 
Epoch JD 2440000.5 Epoch JD 2443000.5 
Mo 244°36047 (1950.0) Mo 59°82706 (1950.0) 
n 0.6058729 w 356°78453 n 0.06871720  w 26°46438 
a 6.420203 Q 32157133 a 5.903273 Q 314. 98038 
e 0.13164292 2 9.46211 e 0.06984021 a 9.12152 
Epoch JD 2441000.5 Epoch JD 2444000. 5 
Mo 302°58791 (1950.0) Mo 123°39207 (1950.0) 
n 0.06123687 w 0°57024 n 0.06747432 w 40°06428 
a 6.374720 2 320.90788 a 5.975544 Q 313.34247 
e 0.13585640 a 9.46080 e 0.06249733 D 9.25203 
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TaBLE ITI. Opposition ephemerides, 1960-1973. 
a (1950.0) 6 p r a (1950.0) 6 p r 
960-61 1962-63 
un 15 2h 47m5 +24° 49’ 6.958 6.185 | Nov 2 62 412 SEBO Op! 
25 DP 54..2 +25 26 12 6 39.5 Sea 09) 6.244 6.921 
ue 5 3 00.4 +26 02 6.758 6.203 22 6 36.7 EG fame 
15 3 06.0 +26 36 Dec 2 6 33.0 ah We) 6.056 6.935 
25 310.9 “LD (OK 6.504 6.221 12 6 28.5 eG Wii 
ug 4 Beto. aia Si 22 6 23.6 +31 29 5.976 6.948 
14 Sais. +28 03 6.218 6.240 | Jan 1 Galons 34 27 
24 3 20.5 98 27 11 @ WGC So 6.021 6.960 
cay. 3 2) SRG) +28 48 5.926 6.258 21 6 09.2 +31 14 
13 SY UA +29 04 31 Gu05e5 +31 04 6.186 6.973 
23 3 20.4 +29 16 5.662 6.276 | Feb 10 6 02.8 Ja) SY 
Yet 3 Sas. +29 22 20 6 01.1 +30 39 6.445 6.985 
13 3 14.8 +1205 23 5.461 6.294 | Mar 2 6 00.6 +30 25 
23 Se HOLY JL Ty 12 6 01.2 +30 12 6.761 6.997 
Nov 2 3 06.0 +29 06 5.357 6.313 22 6 02.9 Be OMEDS 
12 3 01.0 +28 49 Apr 1 6 05.6 +29 45 7.097 7.009 
22 2 56.2 +28 28 5.371 6.331 11 6 09.2 +29 33 
Yec 2 2 51.8 +28 04 21 G ies +29 21 7.417 7.020 
12 2 48.0 +27 39 5.506 6.349 | May 1 6 18.8 +29 09 
22 D452 eT add: 11 6 24.6 SEKI | 7.692 7.032 
fan 1 2 43.5 +26 52 5.743 6.368 Var EPOnS =5 
11 2 42.9 +26 34 oe 
21 De Ry +26 19 6.049 .386 
31 2 45.1 +26 10 es : 
“eb. 10 2 47.9 +26 06 6.387 6.404 | Aug 19 7 33.9 +26 41 7.895 7.084 
20 D516 +26 06 29 7 40.1 +26 27 
Mar 2 256.2 +26 11 6.719 6.422 | Sep 8 7 45.8 =-20) alls 7.685 7.094 
12 S026 DKS Oi 18 7 50.9 +26 04 
22 3007.6 +26 32 7.015 6.440 28 7 55.4 +25 55 7.414 7.104 
Var Shy eds Oct Ss i B.A +25 47 
18 8 01.9 JEN ily 7.105 7.113 
[961-62 28 8 03.8 +25 40 
Jul 20 4 38.0 +30 31 7.222 6.547 | Nov 7 8 04.7 +25 41 6.790 7.122 
30 4 44.4 +30 49 17 8 04.6 +25 44 
Aug 9 4 50.2 +31 07 6.988 6.564 27 8 03.4 +25 50 6.505 7.130 
19 Hae sae 2. +31 24 Dee? 8 01.2 +25 57 
29 4 59.4 +31 40 6.709 6.581 17 i Bysis il +26 06 6.290 7.139 
Sep 8 5 02.6 +31 56 27 7 354.2 +26 15 
18 5 04.8 Se ali 6.410 6.598 | Jan 6 Ue CORE +26 24 6.180 7.147 
28 5 05.8 +32 24 16 7 45.0 +26 32 
Oct 8 See 0527 +32 36 6.124 6.615 26 7 40.3 +26 37 6.194 7.155 
18 5 04.3 +32 46 Feb 5 i SY) +26 39 
28 5 01.9 2132) 53 5.886 6.632 15 7 32.0 +26 39 6.330 7.162 
Nov 7 4 58.4 +32 56 25 ih PSO) +26 36 
17 Aa gies +32 55 5.733 6.649 | Mar 6 7 26.8 +26 30 6.567 7.170 
27 4 49.2 ==82) 50 16 if PSodl +26 22 
Dec 7 4 44.0 +32 40 5.694 6.666 26 (ih PSG: +26 11 6.871 7.177 
17 4 39.0 Fao) 26 Apress) 7 26.6 +25 59 
27 4 34.3 +32 10 5.779 6.682 15 728.0 +25 45 7.201 7.183 
Jan 6 4 30.4 +3452 25 ih hilass +25 30 
edo ea RS) 4-31 33 5.977 6.699 | May 5 7 34.9 +25 14 7,922 7.190 
26 495.4 +31 14 15 7 39.2 +24 56 
Feb 5 4 24.7 +30 57 6.259 6.715 25 7 44.2 +24 36 7.805 7.196 
15 4 25.2 +30 42 Var +2.4 —11 
3 - 4 26.8 +30 29 6.588 6.731 
ar eee o 5 +30 20 
17 4 33.2 30-13 6.925 6.747 | 1964-65 
27 437.8 +30 09 Sep 22 8 57.9 +19 44 7.866 e227 
Apr 6 4 43.1 +30 06 7.238 6.762 | Oct 2 9 02.7 ed Ol N22 
16 4 49:1 +30 05 12 9 07.0 +19 03 7.597 7.231 
26 2155.7 +30 05 7.499 6.778 22 9 10.5 +18 46 
Var +3.1 dike Nov 1 Omiar3 Sg BW 7.284 7.235 
11 Omm5 =18) 22 
1962-63 21 9 16.3 +18 16 6.959 7.239 
Aug 4 6 10.8 +30 24 7.609 6.852 | Dec 1 9 16.4 +18 14 
- 14 617.3 +30 25 11 9 15:5 +18 16 6.658 7.242 
24 6 23.3 ete) 27, 7.389 6.866 21 OF 1357, ai Al 
Sep 3 6 28.6 +30 29 31 9 11.0 +18 30 6.422 7.245 
13 6 33.1 +30 32 7,145 6.880 | Jan 10 9 07.5 +18 41 
23, 6 36.8 +30 36 20 9 03.6 +18 54 6.287 7.248 
Oct .3 6 39.5 +30 41 6.811 6.894 30 S715012 +19 07 
13 ewal.2 +30 47 Feb 9 8 54.8 +19 19 6.275 7.250 
23 6 41.8 +30 54 6.509 6.908 19 8 50.6 +19 30 
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TaBLe III (continued) 
a (1950.0) 6 p r a (Q9S0%0)i26 p r 
1965 1967 
Mar 1 8h 4629 +19° 38’ 6.386 7:252)| Jen: 99 10 526 “£010.56 
11 8 43.8 +19 43 19 10) S5ed +01 43 7.224 7.076 
21 8 41.5 419 45 6.603 7.254 29 10 58.3 +01 26 
31 8 40.1 419 43 Jul 9 11 02.4 +01 04 7.511 7.066 
Apr 10 8 39.6 +19 38 6.892 7.255 19 11 06.4 +00 38 ; 
20 8 40.1 +19: 30 29 14, a2 +00 08 7 Pot 7.055 
30 8 41.5 199 18 7 243 7250 Var +1.9 16 
May 10 8 43.8 +19 04 e) ee a 
20 8 46.8 +18 46 ihe i 
30 8 50.6 +18 26 eat i 
Jun 9 8 54.9 +18 03 7.818 7.258 | Nov 6 12 08.1 —06 38 7.748 6.999 ; 
19 8 59.7 +17 38 16 12 13.0 —07 19 C 
Var aD it Si Dec 6 1225 —08 35 Lh 
16 12 24.8 —09 08 7.184 6.975; 
26 12373 ‘09° 37 
1965-66 Jan 5 12 29.1 10+ 02 6.845 6.962 
Oct 7 10 02.9 +11 44 7.932 7.256 15 12 29.9 —10- 22 ' 
17 10 07.5 46t1 14 25 12 29.9 =10 37 6.511 6.949 © 
Dy 10s Hikes +10 46 7.665 7.254 | Feb 4 12 28.9 =10. 45 4 
Nov 6 10 15.0 +10 21 14 17° 27..0 —10 48 6.233 6.937 
16 10. 17:8 +10 00 7.349 (1232 24 12 24.4 —10 44 
26 10 19.7 +09 43 Mar 5 12: 24.0 =10; 334 6.017 6.923 © 
Dec 6 10 20.9 +09 31 7.016 7.250 15 19-472 —10 19 j 
16 Mie +09 23 25 42 13.2 —10 00 5.922 6.910 
26 10 20.4 +09 20 6.703 7.248 | Apr 4 1209.1 —09 38 4 
Jan.55 10 18.9 +09 22 14 12) 105:3 —09 14 5.948 6.896 
15 10 16.5 +09 29 6.449 7.245 24 12 01.9 08.5 5 
25 10 13.4 +09 40 May 4 1th 35029 08" 30 6.087 6.882 
Feb 4 10 09.7 +09 54 6.291 7.242 14 {157.3 —08 13 7 
14 10 05.7 +10 10 24 11. 56.2 07) 58 6.314 6.868 
24 10 01.6 +10 27 6.254 7.239 | Jun 3 1 558 —07 50 L 
Mar 6 9 57.6 410 43 13 11 56.6 —07 47 6.596 6.854 
16 9 53.9 +10 58 6.342 7.236 23 ese 3 —07 50 h 
26 9 50.8 eit 14 jal 33 12 00.5 ieee 6.898 6.840 
Apr 5 9 48.5 4-11 20 6.537 7.232 13 12 03.6 —08 10 ; 
15 9 46.9 4105-95 23 12 07.4 —08 29 7.186 6.825 
25 9 46.2 sediee2s 6.809 7.228 | Aug 2 12 -t,.8 —08 52 
May 5 9 46.5 edhe ge 12 t2> A6.7 —09 19 7.435 6.810 © 
15 9 47.6 eS nay! el 7.223 Var 04 —18 
25 9 49.5 +11 02 i 
Jun 4 O SP +10 46 7.435 7.219 Si 
14 9 55.6 +10 26 ape | 
24 0) 505 +10 02 7799 7.214 | Dec 10 163 PADS —17 24 7.292 6.717 j 
july 10 04.0 +09 35 20 13 34.1 —18 03 q 
14 10 08.9 +09 06 7.956 7.208 30 13 38.2 =i8 40 6.994 6.701 | 
Var aig —16 Jan 9 13, 41567 Fo 219.43 
19 13 44.2 —19 43 6.659 6.684 © 
1966-67 29 13 45.9 —20 09 : 
Feb 8 13 46.6 —20 29 6.320 6.668. 
Oct 22 11 04.7 +02 48 7.889 147.178 18 13 46.3 —20 45 
Nov 1 11 09.4 +02. 11 28 13 44.9 +20 54 6.014 6.651 
11 ft eis +01 37 7.626 7.171 | Mar 10 13 42.7 =20° 57 q 
21 11 470 +01 06 20 13 39.6 —20 54 OM i 6.634 
Dec 1 11 19.9 +00 38 7.310 7.163 30 13 35.8 —20 44 gq 
11 (e220 +00 13 Apr 9 13, 31.6 —20 29 5.642 6.617 
21 11 23.4 —00 04 6.972 7.156 19 13273 —20 09 
31 it $2308 200,48 29 13) 23.0 —19 46 5.624 6.600 
Jan 10 11 23.4 —00 26 6.648 7.148 | May 9 13 19.1 S19, 23 
20 Mpa —00 29 19 13° 15.8 S57 522 6.583 
30 11 19.9 —00 26 6.377 7.140 29 13 13.3 —18 34 
Feb 9 WR ES =005 17 Jun 8 ES ORTINE =18° 44 5.916 6.566 
19 11 13.6 —00 04 6.197 7.132 18 13 11.0 =17. 58 
Mar 1 11 09.7 +00 13 28 13 SA 3 =—17 47 6.176 6.548 
11 11 05.7 +00. 32 6.133 7,423) yale os 13 12.5 il fn AP. 
21 At O1eS +00 53 18 rae ey Aare) 6.468 6.531 
31 10 58.1 +01 13 6.193 7.114 28 13.7 47..8 =i7 48 
Apr 10 10 55.0 +01 31 Aug 7 13 21.6 —17_ 59 6.759 6.513 
20 10 °52"5 +01 46 6.364 7.105 17 13 26) 1 S18 245 
30 10 50.8 +01 58 27 iseas 3 18.36 7.019 6.495 
May 10 10 50.0 +02 05 6.616 7.096} Sep 6 1337.1 —19 00 | 
20 10 50.0 +02 07 16 13 43.3 219° 28 7.22% 6.477 
30 10 50.9 +02 04 6.915 7.086 Var +2.5 —18 


COMET SCHWASSMANN-WACHMANN Dit 


TABLE III (continued) 


a (1950.0) 6 p r a (1950.0) 6 p r 
969-70 1972 
Yec 25 14h 55™2 —25° 33’ 7.050 6.386 | Feb 13 18h 57™3 —28° 45’ 6.421 5.711 
an ° 4 15 01.4 —26 08 23 19 05.5 —28 35 
14 ily OAaal —26 Al 6.774 6.368 | Mar 4 1ADTAS.0 —28 25 6.171 5.697 
24 LS aan —27 13 14 19 19.8 —28 15 
‘eh «3 TS LO. 2 —27 43 6.453 6.350 24 19 25.8 —28 08 5.871 5.683 
13 15 19.3 —28 09 Apr 3 19 30.8 —28 03 
23 15 21.4 —28 33 6.114 6.331 13 19 34.7 —27 59 5.547 5.670 
arn 5° 5m 224 —28 33 23 1993755 —27 57 
15 IGS BRS) —29 10 5.791 6.313 | May 3 19 39.0 —27 59 5.227 5.657 
25 S20 —29 21 13 19 39.3 —28 03 
\pr 4 15 18.6 —29 27 5.521 6.294 23 19 38.3 —28 08 4.944 5.645 
14 sho Cassel —29 27 Apeney 2 19 36.0 —28 15 
24 1 edc0 —29 21 Sos 6.276 12 19 32.6 —28 22 4.732 5.633 
May 4 15 06.4 —29 10 ; D2. 19 28.3 —28 28 
14 15 01.6 —28 53 5.263 OF2578 |e jul) 92 19002373, —28 33 4.618 5.621 
24 14 56.8 —28 32 12 19 18.0 —28 36 
qi) 3. 14 52.6 —28 09 5.305 6.238 22 1927S —28 35 4.619 5.609 
13 14 48.9 —27 45 Aug 1 19 07.9 —28 31 
Aes 14 46.3 —27 21 5.454 6.220 11 19 03.8 —28 23 4.733 5.598 
mle 3 14. 44.5 —26 59 Mii 19 00.7 —28 13 
13 14 44.0 —26 41 5.683 6.201 31 18 58.7 —28 ° 00 4.941 5.587 
23 14 44.5 —26 27 Sep 10 18 58.0 —27 45 
lug 2 14 46.2 —26 17 5.960 6.183 20 18 58.7 —27 28 5.214 Sel 
12 14 48.9 269 ale) 30 19 00.7 —27 09 
22 14 52.6 —26 12 6.252 6.164 | Oct 10 19 03.8 —26 50 5.516 5.567 
ep) iil yf ae —26 16 20 19 08.1 —26 30 
11 iS: OPS —26 24 6.529 6.146 30 19 13.4 —26 08 5.816 5.558 
21 15 08.8 —26 36 Nov 9 19 19.6 —25 44 
Met 1 15 95;,.6 —26 51 6.765 6.127 19 19 26.6 —25 19 6.084 5.548 
Var +3.3 —13 Var +4.2 +15 
1971 1973 
an 29 16 55.9 —30 56 6.568 6.018 | Mar 19 21 09.7 —18 19 6.180 5.502 
‘eh «8 17 02.8 —31 14 29 21 16.9 —17 42 
18 17 08.9 —31 31 6.273 6.001 | Apr 8 21 23:4 —17 07 5.930. 5.496 
28 Lip: —31 47 18 21 29.3 —16 34 
Mar 10 Wie tees —32 03 5.945 5.983 28 WL Shee —16 04 5.635 5.490 
20 sy DHS —32 18 May 8 2A 38e5 —15 38 
30 M2304 —32 32 5.613 5.965 18 21 41.6 —15 16 5.320 5.484 
\pr 9 Tee 2420 —32 46 28 PARAS ee —14 59 
19 Li 32323 —32 57 Sole 5.948 | Jun 7 21 44.6 —14 47 5.013 5.479 
29 1721.53 —33 06 17 Dil LYS —14 40 
May 9 Ye wilciaa —33 12 5.075 5.931 27 Di e43e1 —14 39 4.745 5.475 
19 17 14.0 —33 14 athe 7 21 40.6 —14 42 
29 fe 509-1 —33 11 4.933 5.914 17 DA e372 —14 50 4.551 5.471 
nn «68 ive NORD —33 04 27 Doon —15 O1 
18 16 58.6 —32 52 4.904 5.897 | Aug 6 21 28.3 —15 13 4.458 5.467 
28 106s 53.0 —32 37 16 Pi GR) —15 25 
ul 8 16 49.3 —32 18 4.988 5.880 26 PA MESES) 59 Si 4.482 5.464 
18 16 45.9 —31 58 Sep, —5 DA at, —15 47 
28 16 43.6 —3t 38 Kai Wal 5.863 15 21 10.5 —15 53 4.619 5.461 
Lug 7 164256 —31 19 25 21 07.8 —15 55 
17 16- 42.8 —31 Ol 5.423 S284 Oct 5 21 06.3 —15 54 4.849 5.458 
27 16 44.3 —30 45 AUS; 21 06.0 —15 48 
ep 66 16 47.1 —30 32 5.713 5.831 25 21 06.9 —15 36 5.139 5.456 
16 16 51.0 —30 21 Nov 4 21 08.9 —15 20 
26 16 55.9 —30 12 6.006 5.815 14 Pak ail —14 59 5.455 5.455 
pet. 16 17 01.8 —30 05 24 2A e16e3 —14 34 
16 17 08.6 —30 00 6.275 5.799 | Dec 4 21 2s —14 04 5.763 5.454 
26 ge elo.t —29 56 14 Dil BIE —13 30 
Nov 5 Hime. 2 —29 352 6.495 5.784 24 PAL Bi Saeh —12 53 6.035 5.453 
Var +4.0 —3 Var +3.7 +28 


The variation given at the end of each opposition 
sphemeris is the variation-per-ten-days for each coor- 
linate in the neighborhood of the opposition. The 
esiduals in the left-hand columns of Table I represent 
he deviations from the preliminary orbit. The right- 
1and columns give the representation of each obser- 


vation from the adopted solution. Each residual enclosed 
in parentheses ( ) was excluded from the final least- 
squares solution. 


REFERENCE 
Herget, P. 1947, Astron. J. 53, 16. 
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Measurement and reduction of plates taken with the Sproul refractor over the interval 1920-1959 yield 
+”038+"005 (p.e.) for the relative parallax and +"008+"010 (p.e.) for the semi-axis major of the photo- 
centric orbit of 42 Comae. The relative mass of the B component is .507+.014 (p.e.), whereas the individual 
masses are poorly determined due to the small value of the parallax: Ma=1.6 ©. Mg=1.7 © each witha 


p.e. of +.45 ©. 


HE visual binary 42 Comae=ADS 8804, 135.1", 
+18°3’ (1900) has a period 25%83 and semi-axis 
major "672 (Pavel 1949). The combined visual mag- 
nitude is 4.30 and the Am,,=.02 according to Eggen 
(1956). The spectrum is F4, Mt. Wilson classification ; 
the radial velocity is —17.7 km/sec (Wilson 1953). The 
inclination of the orbit is close to 90°, the greatest 
separation never exceeds 0"6, so that the image of AB 
always appears circular on the Sproul photographs. 
The procedure described in previous Sproul publi- 
cations has been followed (van de Kamp and Lippincott 
1949). A total of 999 exposures on 300 plates covering 
98 nights with a total weight 208 were measured by 
Mrs. Mary Jackson on the St. Clair-Kasten measuring 
machine. The brightness of 42 Comae was reduced 6.7 
mag. by means of a rotating sector to equal the average 
magnitude of the reference system of three stars. 
Table I gives information about the reference stars 
as well as the dependences for 1940 and their annual 
changes. The coordinates of the standard frame and of 
42 Comae hold for the approximate epoch 1940. 
The nightly means were represented by the equations 
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Separate solutions in right ascension and declination 
yield the following results (the scale factor is 1 mm’ 
eS ISY EWA) 
mm pres wt 
Cz= — 3.59069 
Cy = +4.29875 
Uz= —0.0220270= —0%415649-+07000808 5894 
by = +0,0085943 = +0.162174 +0.000530 5971 


m-= +0.002131 =+0.04021 +0.00723 73.62. 

m= +0.002322 =+0.04383 +0.00804 26.12 

a.=—0.01691 =—0.3191 -£0.0574 1.16 

ay=+0.00121 =+0.0228 -+0.0077 26.34 
pe t= +00033  =40.062 | 
p.e. 1y=+0.0022 =+0.041 


A least-squares solution for + and a was made com- 
bining the material in right ascension and declination 
with the following results: 


of conditions f: 35 D i “i 
X=crtudttaPet an aged 
Veen Cy= + 4.29876 
2= —0"741903+-0"00047 12202 
where the symbols have the usual meaning; ¢ is counted f 
from 1940.00. Hy=+0.16286 +0.00066 6126 
The orbital factors were computed from Pavel’s a=+0.0082 +0.0097 27.55 
(1949) orbit: m= +0.0375 +0.0054 90.54 
P=25383, 1= 89°94, p.e. 1=+0.0513 
= 494 = DNS . 
§ i et While the extent of the reference configuration is 
T=1911.61, Q=11°93. large in right ascension, it does not fully account for 
TABLE I. Reference stars. 
: Dep. 
No. Diameter Mpy Spectrum Wa Vs 1940 AD/year 
1 "7093 11.3 GO 76,31 21.73 305 00025 
2 .104 ileal KO +25.57 —15.53 . 166 + ooo 
3 .099 TE F8 +50.80 +17.26 -439 + .00016 
42 Com .110 (10.98) F4 — 3.59 + 4.30 


| 
| 


42 COMAE 


TABLE II. Normal points of residuals. 


Number Total Vz, 2, 


Epoch of nights weight unit .0001 mm 

1921.22 5 3 +7 —10 
22.80 6 34 —36 + 8 
26.67 4 3 —67 +25 
38.91 4 5 —2 — 3 
40.64 4 7 +10 — 3 
43.91 4 8 +43 +5 
46.85 9 193 +15 0 
48.97 8 21 +16 —16 
51.01 5 12 — 6 —10 
52.25 5 i, —21 are 
Saya 5 13 + 2 "+13 
54.21 7 ilgf — 1 — 6 
55.87 12 32 0 0 
57223 6 13 —27 +16 
58.24 9 27 0 0 
59.18 5 12 —7 — 3 


_ the large value of the probable error of unit weight 


given by the least-squares solution (Lippincott 1957). 
Table II gives normal points of the residuals found from 
the combined solution. These residuals also show greater 
scatter in right ascension than in declination. 

This solution for parallax supersedes the early Sproul 
determination (Pitman 1928), which was based on part 
of the present material. The reduction to absolute is 
made in the usual manner (Lippincott 1958) and leads 
tO Tabs =+0”042+0"005. Combined with the values 
for the other published determinations (Jenkins 1952) 
the weighted mean value for the absolute parallax is 
+0"051+0"004. With this value, the semi-axis major 
a becomes 13+1 a.u., yielding Ma+Mp=3.340.8 ©, 
the p.e. in both cases arising solely from the p.e. of the 
parallax. Adopting Eggen’s value Am=0.02, B=0.495; 
we derive the mass ratio 
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B= (a/a)+8=+0.507+0.014 (p.e.) 


taking into account only the p.e. of a (van de Kamp 
1945). 

Thus we find Ma=1.6 ©, Mg=1.7 © each with a 
p-e.=+0.45 © due to the large percentage error in the 
parallax. The absolute visual magnitudes of the two 
components are 3.58 and 3.60, the absolute bolometric 
magnitudes 3.53 and 3.55. Considering the probable 
errors of the masses, the two components do not deviate 
from the mass-luminosity relation. 

Pavel’s orbit gives i=89°56’+6’. For eclipses to take 
place the inclination would have to be within 1’ of 90°. 
Assuming for the moment this to be the case, the next 
eclipse would occur very close to periastron lasting 
roughly two hours; the uncertainty of the parameters 
precludes a more accurate timing of smallest angular 
separation so that without continuous 24-hour obser- 
vations for several days the occurrence or absence of an 
eclipse could not be established. 
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NOTICES 


Observing Time on Kitt Peak ; 


HE Kitt Peak National Observatory was estab- 
lished to provide American astronomers and 
qualified graduate students with modern telescopic and 
instrumental equipment for astronomical research. Also, 
the Board of Directors of AURA (Association of Uni- 
versities for Research in Astronomy) has set the policy 
that about 60% of the observing time on Kitt Peak 
shall be made available to visiting astronomers and 
graduate students. 

No charge is made for telescope time, but visiting 
astronomers will be billed for food and lodging while 
working on Kitt Peak. 

At the present time, two moderate-sized, Cassegrain 
focus telescopes are in operation on Kitt Peak: a 16-inch 
and a 36-inch. Auxiliary equipment consists of photo- 
electric photometers for both telescopes, and a grating 
spectrograph for the 36-inch. 

Since a number of inquiries and requests for observing 
time have already been received, an Observers’ Infor- 
mation Leaflet and an Observing Time Request Form have 
been prepared for those persons who are interested in 
using the Kitt Peak facilities. Copies may be obtained 
by writing to: Office of the Director, Kitt Peak National 
Observatory, 950 North Cherry Avenue, Tucson, 
Arizona. 

It is AURA policy that observing-time requests be 
granted on the basis of scientific merit of the proposed 
program. In most cases, the Observatory staff reviews 
proposals and makes recommendations regarding as- 
signments; occasionally, the advice of outside referees 
may be sought. Visiting astronomers are expected to do 
their own observing on Kitt Peak, although they may 
bring a qualified assistant if their program réquires it. 

Subject to existing Federal regulations, qualified 
foreign astronomers may use the Observatory’s facilities. 


; 
f 


In cases of U.S. and foreign programs judged to be of 
equal merit, however, preference will be given te 
astronomers resident in the U.S., because of the national 
character of the Observatory. 

To the extent that funds are available, the Obsex 
vatory’s graduate student program pays for travel 
and subsistence expenses of qualified, advanced-degree 
U.S. students. 

Although the resident staff wishes to be as helpful as 
possible to visitors, its responsibility is limited to giving 
advice and initial assistance for the efficient use of the 
instruments. The Observatory staff and assistants will 
not do observing for visitors except under emergency 
conditions. 

In addition to the telescopes on Kitt Peak, there are 
at the Tucson headquarters a number of laboratory 
instruments, such as an electronic digital computer, a 
Hilger microphotometer, a Cary monochromator and 
spectrometer, measuring machines for direct photo- 
graphs and spectrograms, a spectrocomparator, and 
binocular microscope viewers. The main Observatory 
library, electronics and photographic laboratories, and 
the instrument and optical shops are also located in the 
Tucson headquarters buildings. 

Depending upon availability of office space, astron- 
omers on sabbatical or other long-term leave may work 
at the Observatory while in residence in Tucson. 

The Observatory staff welcomes inquiries from those 
astronomers and graduate students who have research 
programs that may merit the use of the Kitt Peak 
National Observatory facilities. 


N. U. MAyatyi 

Observatory Director 

Kitt Peak National Observatory, 
Tucson, Arizona 
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The Council and Members of the American Astro- 
‘ nomical Society have studied the proposed project, 
| West Ford, as described in a series of papers in the 
_ Astronomical Journal (Astron. J. 66, 105-118, 1961). If 
| the project is successfully carried out according to the 
published plans, it will result in a belt of orbiting di- 
poles around the earth which will be barely detectable 
| by astronomers and radio astronomers with present 
techniques when observations are carried out from the 
ground. The experimental belt could, however, interfere 
seriously with delicate and important observations from 
satellite observatories scheduled to begin operation 
within the next two to three years, unless the lifetime 
of the belt is sufficiently limited. Consideration also of 
the rapid gains that can confidently be expected to 
occur in the sensitivity of astronomical measurements 
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Resolution Adopted by American Astronomical Society, June 20, 1961, 
at Nantucket, Massachusetts 


during the next few years makes it imperative to guard 
against the irreversible contamination of space. The 
Society is even more gravely concerned at the harm that 
a future operational belt or belts might cause to astro- 
nomical observations, and also calls attention to the 
interference that can be expected from proposed systems 
of passive and active communication satellites. 

In the light of the foregoing considerations, the So- 
ciety resolves as follows: 


1. The American Astronomical Society strongly op poses 
any contamination of space that is detrimental to the 
conduct of basic scientific research. 

2. The Society reaffirms its previous resolution in sup- 
port of world-wide frequency allocations for radio as- 
tronomy as a matter of great and increased urgency. 


ERRATA 


Astron. J. 53, 6. 
For Mo=350°56422 read Mo=353°56422. 


(Communicated by the author, Dr. Paul Herget.) 


Astron. J. 65, 385. 


At JDE 2434980.5: for Mo=1°41280 read 1°841064. 
At JDE 2437420.5: for My>=0°81742 read 2°32024. 
No other parts of the work are affected. 


(Communicated to the author by Dr. J. G. Porter.) 


Astron. J. 63, 338. 
Table 1. No. 1: for +1°990 read +0°990. 
No. 2: for +1°989 read +0°989. 
Delete stars Nos. 23, 26, and 28. 
No. 34 for +5'41™6 read +5525™6. 


(The authors express their thanks to Mr. Henry Giclas of the Lowell Observatory for calling their attention 


to these errors.) 


Astron. J. 66, 46. 


‘The final paragraph of the abstract by Hindman, McGee, Carter, and Kerr, reading ‘The peak . . . per 
year,” has been misplaced through editorial accident. This paragraph should conclude the abstract by McGee 


and Murray on page 48. 


i 
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NEW BOOKS RECEIVED 


Mécanique Céleste, by F. Tisserand. Published by Gauthier-Villars, Paris, 1960. Vol. I, 50 N.F.; Vol. II, 60 N.1 


Photographically reproduced in paperback binding (63X93 inches) from the first two volumes which fir: 
appeared in 1889 and 1891, respectively. According to the Preface of A. Danjon, all known typographical erro: 
have been corrected ; otherwise the volumes have not been revised. 

Tisserand’s Mécanique Céleste is deservedly famous for its lucid exposition. Even after 70 years, there is muc 
in these two volumes on which it would be difficult to improve. 

Volumes III and IV are still available in the original edition. 


Methods of Celestial Mechanics, by Dirk Brouwer and Gerald M. Clemence. Published by Academic Press Inc 
New York and London, 1961. Price $15.50. 


Workers in the field have long felt the need for an authoritative, modern treatise which is simultaneously of us 
to advanced students and research workers. The present volume fills this need most successfully. Further tk 
treatment is built up sufficiently from first principles to permit use of the book for self-study, as an introductor 
text for graduate instruction, or as a reference volume. The book establishes a solid balance between theoretic: 
treatment and practical applications of the subject. Particularly helpful are worked-out numerical examples 1 
the chapter on numerical integration of orbits.* 


An Introduction to Celestial Mechanics, by T. E. Sterne, Series No. 9, Interscience Tracts on Physics and Astronomy 
Published by Interscience Publishers, Inc., New York. Price $4.50 (cloth) and $2.50 (paper). 


A competently written introduction to the basic topics of celestial mechanics, followed by applications to tl 
problem of artificial satellite motion (including effects of atmospheric drag), and a chapter on numerical integr 
tion. The author has succeeded in presenting this material clearly and attractively. Yet many readers may fe 
the need to consult more elaborate treatises. This applies, for example, to the indication of the character of gener 
planetary and lunar theories (pp. 108-115), but also to some sections in the earlier chapters. This is, of cours 
unavoidable in such a small book and does not detract from its considerable merits. 


\ 


* The senior Astronomical Journal editor wishes to disassociate himself from responsibility for the above notice which was kind 
contributed by another reviewer. 
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